F/0  20/1 


AD-A116  309 
UNCLASSIFIED 


u  I 


TEXAS  UNIV  AT  AUSTIN  APPLIED  RESEARCH  LABS 
ACOUSTIC  REFLECTION  FRON  ANISOTROPIC  SEDIMENTS. CUI 


JUN  S2  DM  OAKLEY 
ARL-TR-Oi-20 


N0001A-00-C-0990 
ML 


•  FILE  COPY 


j.v*; 


•^s-yg-  yyr  VIT*'****  **T*-t"J 


Mti'.'.wwflfinwii^fy  upp! 


*\ 


O 

00 

CO 


.'A  Vi-  V-V* 


■■  V>^] 


.-a 


ARL-TR-82-28 


Copy  No.  3  / 


ACOUSTIC  REFLECTION  FROM  ANISOTROPIC  SEDIMENTS 


iH 

< 


David  W.  Oakley 


APPLIED  RESEARCH  LABORATORIES 

THE  UNIVERSITY  OF  TEXAS  AT  AUSTIN 
POST  OFFICE  BOX  S029,  AUSTIN.  TEXAS  7S712-S029 


1  June  1982 
Technical  Report 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


Prepared  for: 

OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VA  22217 


82  07  02  027 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PACE  flWiMi  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

L  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

4t~  A  //6  2 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtltla) 

ACOUSTIC  REFLECTION  FROM  ANISOTROPIC  SEDIMENTS 

s.  type  of  report  ft  period  covered 

technical  report 

6  PERFORMING  ORG.  REPORT  NUMBER 

ARL-TR-82-28 

7.  AUTHORf.J  . 

Davi d  W.  Oakley 

8.  CONTRACT  OR  GRANT  NUMBERS) 

N00014-80-C-0490 

9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Applied  Research  Laboratories 

The  University  of  Texas  at  Austin 

Austin,  Texas  78712-8029 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 

AREA  A  WORK  UNIT  NUMBERS 

Item  0015 

It.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Office  of  Naval  Research 

Department  of  the  Navy 

Arlington,  VA  22217 

12.  REPORT  DATE 

1  June  1982 

13-  NUMBER  OF  PAGES 

14.  MONITORING  AGENCY  NAME  ft  ADDRESS^)/  dlllarant  Iron l  Controlling  Ollier) 

IS.  SECURITY  CLASS,  (ol  thla  taport) 

UNCLASSIFIED 

15«.  DECLASSIFICATION'  DOWNGRADING 
SCHEDULE 

l«.  DISTRIBUTION  STATEMENT  (ol  thla  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  ebatrect  entered  In  Block  20,  ft  different  from  Report ) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  reveree  aide  if  neceaamry  end  Identify  by  block  number) 

2CL?  ABSTRACT  (Continue  on  reveree  aide  If  neceaaery  end  Identity  by  block  number) 

Measurements  of  an  excess  of  horizontal  to  vertical  velocity  (averaging  10" 
for  shales)  have  been  made  in  marine  sediments  and  in  sedimentary  rocks.  In 
this  paper,  the  proposed  causes  of  such  anisotropies  are  reviewed,  and  the 
reflection  coefficient  for  incidence  from  water  onto  a  solid  sea  floor  is 
derived  for  plane  wave  solutions  of  the  anisotropic  wave  equations.  The 
possibility  of  backward  traveling  fronts  of  pressure  amplitude  at  low  grazing 
angles  is  noted.  A  propagation  matrix  formation  is  established  to  generalize 
the  reflection  coefficient  for  the  case  of  multiple  layers.  Bottom  loss^ 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWi»n  Data  Entarad) 


DD  I  JAN  73  1473  EDITION  OF  I  NOV  «5  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE (Whmti  Dmtm  Entmrmd) 


20.  (Cont'd) 

calculations  are  made  which  suggest  that  anisotropy  can  increase  the 
reflectivity  of  the  ocean  floor. 


COPY 

.  INSPECTED  . 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  P AGEfWftmn  Dm tm  Entmrmd) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES  vii 

LIST  OF  FIGURES  viii 

CHAPTER  I  INTRODUCTION  1 

CHAPTER  II  OCCURRENCE  AND  CAUSES  OF  VELOCITY  ANISOTROPY  3 

A.  Independent  Elastic  Parameters  for  Transverse  11 

Isotropy 

B.  Equations  of  Motion  for  Transversely  Isotropic  21 

Media 

C.  Plane  Wave  Solutions  to  the  Equations  of  Motion  25 

D.  Reflection  of  a  Plane  Wave  from  an  Anisotropic  29 

Solid 

E.  Reflection  from  a  Layered  Anisotropic  Solid  32 

CHAPTER  IV  PROPERTIES  OF  THE  REFLECTION  COEFFICIENTS  41 

A.  Effects  of  Anisotropy  on  Wave  Numbers  45 

B.  Effects  of  Anisotropy  on  Fluid-Solid  Interface  47 

Reflection  Coefficient 

C.  Parameter  Studies  on  the  Constant  C^  56 

D.  Effects  of  Anisotropy  on  Reflection  from  a  Typical  62 
Marine  Sediment 

CHAPTER  V  CONCLUSIONS  71 

BIBLIOGRAPHY  75 

iii 


LIST  OF  TABLcS 


E 


fci 

H 

f  ■ 

i 

r  i 
f  i 


Table 

I 

II 


III 


IV 


Representative  Velocity  Anisotropies  from  DSDP  Data 
Anisotropic  Elastic  Parameters  for  Green  River  Shale, 
Measured  under  8000  psi  Confining  Pressure 
Anisotropic  Elastic  Parameters  for  Green  River  Shale, 
Measured  under  1000  psi  Confining  Pressure 
Anisotropic  Elastic  Parameters  for  Williston  Basin 
Cretaceous  Shale 

Possible  Interpretations  of  Calculated  Compressional 
Velocity  for  Anisotropic  and  Isotropic  Sediments 


Page 

5 

42 

43 

44 

65 


LIST  OF  FIGURES 


Figure  Page 

1  The  Vertical  Wave  Numbers  in  the  Complex  Plane  46 

2  Bottom  Loss  for  Fast  Green  River  Shale  49 

3  Phase  of  Reflection  Coefficient  for  Fast  Green  50 

River  Shale 

4  Bottom  Loss  for  Slow  Green  River  Shale  51 

5  Phase  of  Reflection  Coefficient  for  Slow  Green  52 

River  Shale 

6  Bottom  Loss  for  Williston  Basin  Shale  53 

7  Phase  of  Reflection  Coefficient  for  Williston  54 

Basin  Shale 

8  Gradual  Reducti on  for  Beta  to  1 ;  Effects  on  Bottom  59 

Loss  for  Fast  Green  River  Shale 

9  Gradual  Reduction  of  Beta  to  1;  Effects  on  Bottom  60 

Loss  for  Slow  Green  River  Shale 

10  Gradual  Reduction  of  Beta  to  1;  Effects  on  Bottom  61 

Loss  for  Williston  Basin  Shale 

11  Gradual  Reduction  of  Beta  to  1;  Effects  on  Phase  of  63 
Reflection  Coefficient  for  Williston  Basin  Shale 

12  Bottom  Loss  from  Soft  Silt-Clay  Sediment;  Beta  =  1  67 

13  Bottom  Loss  from  Soft  Silt-Clay  Sediment;  Beta  =  0.75  68 

14  Bottom  Loss  from  Soft  Silt-Clay  Sediment;  Beta  =0.5  70 


vii 


CHAPTER  I 


INTRODUCTION 

Customary  studies  of  acoustic  bottom  interation  model  the 
sea  floor  as  a  series  of  flat,  isotropic  layers  underlaid  by  a  flat,  semi¬ 
infinite  homogeneous  isotropic  substrate.  However,  it  has  always  been 
recognized  that  this  model  represents  an  extreme  idealization  of  the  actual 
bottom  environment. 

The  assumption  of  isotropy  of  sediments  has  always  been  known  to 
be  incorrect  on  a  microscopic  scale.  The  mineral  grains  composing  the 
various  sediments  are  inherently  anisotropic.  Acousticians  have  assumed 
that  because  these  anisotropic  grains  are  randomly  distributed,  the  sedi¬ 
ment  therefore  presents  overall  isotropy. 

Unfortunately,  experimental  evidence  does  not  bear  out  this 
assumption.  Petroleum  geophysicists  have  reported  velocity  anisotropies 
from  their  well  log  data  for  many  years  and  the  work  of  the  Deep  Sea 
Drilling  Project  (DSDP)  has  confirmed  that  velocity  anisotropies  do, 
indeed,  exist  in  marine  sediments. 

It  is  the  purpose  of  this  paper,  therefore,  to  examine  the  causes 
of  velocity  anisotropies  in  sediments  and  the  properties  of  those  sedi¬ 
ments,  and  to  develop  the  mathematical  theory  of  reflection  from  such  sedi¬ 
ments.  The  particular  mathematical  framework  which  has  been  chosen  is  that 
for  transversely  isotropic  layers,  that  is,  layers  which  have  a  single 
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unique  axis,  and  rotational  symmetry  about  that  axis.  Propagation  in 
such  a  medium  will  be  examined,  and  a  parameterized  version  of  a  reflec¬ 
tion  coefficient  model  developed.  The  theory  utilized  in  the  work  is 
based  on  that  of  Stone! ey.  His  derivation  of  the  equation  of  motion  is 
followed  and  used  as  a  basis  for  our  derivation  of  the  reflection 
coef f i ci ent . 

Chapter  II  of  this  work  is  a  discussion  of  observed  anisotropies 
in  marine  sediments  and  the  mechanisms  considered  responsible  for  such 
anisotropies.  The  theoretical  form  of  the  reflection  coefficient  is 
derived  in  Chapter  III,  and  in  Chapter  IV  the  results  of  numerical  cal¬ 
culations  in  the  formulation  are  compared  with  those  of  standard  isotropic 
theory. 
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CHAPTER  II 


OCCURRENCE  AND  CAUSES  OF  VELOCITY  ANISOTROPY 

The  most  commonly  noticed  anisotropy  in  marine  sediments  is  an 
excess  of  the  parallel  shear  and  congressional  velocities  over  the  veloci¬ 
ties  normal  to  the  plane  of  bedding.  Hamilton  (1970)  concludes  from  data 
obtained  by  LeRoy  (1950)  and  Uhrig  and  Van  Melle  (1955)  that  shale 
exhibits  an  average  anisotropy  of  horizontal  over  vertical  velocities 
of  about  10%,  though  the  excess  may  run  as  high  as  40%. 

Milholland,  Manghnani ,  Schlanger,  and  Sutton  (1980)  state  a 
measure  for  this  anisotropy, 


V 


where  subscripts  h  and  v  refer  to  the  directions  parallel  and  normal, 
respectively,  to  the  plane  of  bedding,  V  refers  to  velocity,  and  super¬ 
scripted  bars  refer  to  averages.  This  same  formula  produces  shear  and 
congressional  anisotropies,  A  and  A  ,  but  it  must  be  remembered  that  the 
propagation  of  a  shear  wave  along  a  given  direction  is  also  dependent  up¬ 
on  its  polarization.  Six  measurements  must  therefore  be  taken  to  char¬ 
acterize  V  for  shear  waves,  measurements  along  each  axis  being  taken  with 
polarizations  both  parallel  and  perpendicular  to  the  plane  of  bedding. 

It  is  this  definition  of  anisotropy  which  is  used  throughout  this  paper. 


Mi  1  hoi  1  and ,  Manghnani ,  Schlanger,  and  Sutton  (1980)  also  define 


a  second  parameter, 

Vhh  .  w+w 

A  _  s  s 


which  measures  the  anisotropy  of  horizontally  propagating  waves  (first 
superscript)  with  polarization  (second  superscript)  parallel  and  normal 
to  the  plane  of  bedding.  This  quantity  is  of  more  value  than  A$  from  a 
theoretical  viewpoint,  as  it  can  be  related  to  the  elastic  constants  of 
the  sediment. 

In  oozes  and  chalks,  Manghnani,  Schlanger,  and  Milholland  (1980) 
report  variations  in  compressional  anisotropy  from  -4%  to  8%,  and  varia¬ 
tions  in  shear  anisotropy  from  -8%  to  6%.  In  carbonate  limestones, 
siliceous  limestones,  and  chert  below  1020  m  they  have  determined  the 
compressional  anisotropy  to  be  8%  ±4%,  and  the  shear  anisotropy  to  be 
6%  ±4%. 

In  a  later  report  from  this  same  group  of  geophysicists 
(Milholland  et  al . ,  1980),  data  from  the  DSDP  is  used  to  obtain  a  more 
detailed  picture  of  the  velocity  anisotropies  of  marine  sediments.  These 
results  are  reproduced  in  Table  I. 

It  is  to  be  noted  from  Table  I  that  carbonate  limestones, 
siliceous  limestones,  and  cherts  have  similar  velocity  anisotropies.  This 
similarity  would  seem  to  imply  that  silica  content  does  not  contribute 
significantly  to  velocity  anisotropy.  On  the  other  hand,  clay  enrichment 
of  calcareous  sediments  increases  velocity  anisotropies,  whereas  non- 
calcareous  clays  show  a  much  smaller  degree  of  anisotropy.  Velocity 
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TABLE  I 

REPRESENTATIVE  VELOCITY  ANISOTROPIES 


FROM  DSDP  DATA 

A_ 

A, 

A 

P 

s 

so 

Sites  288,  299 

ill 

ill 

ill 

Carbonate  Ooze 

2.0 

Carbonate  Chalk 

2.9  (3.3) 

0.5  (4.7) 

2.7  (4.2) 

Carbonate  Limestone 

8.0  (4.7) 

6.2  (4.9) 

7.9  (4.7) 

Siliceous  Limestone 

8.6  (3.7) 

5.6  (3.4) 

5.9  (4.0) 

Chert 

8.3  (2.3) 

6.3  (2.0) 

6.2  (3.7) 

Basalt 

2.5  (1.0) 

A_ 

A 

A 

P 

s 

so 

Site  210 

ill 

ill 

ill 

Cl  ay- bearing  Ooze  &  Chalk 

6.2  (1.4) 

6.63  (3.7) 

9.6  (4.5) 

Clay-rich  Chalk 

6.8  (2.6) 

14.2  (5.9) 

9.2  (6.4) 

Calcic-rich  Chalk 

8.1  (1.2) 

17.4  (5.0 

23.0  (5.6) 

Non-calcareous  Clay 

3.3  (1.2) 

-0.4  (2.4) 

12.2  (6.7) 

The  numbers  in  parentheses  represent  the  standard  deviation  for  the 
corresponding  sediment  type. 


anisotropies  in  clay-bearing  to  clay-rich  oozes  and  chalk  are 
approximately  twice  as  large  as  in  calcareous  oozes  and  chalks.  Further¬ 
more,  shear  anisotropies  and  shear  orientation  anisotropies  tend  to  be 
greater  than  compressional  anisotropies  for  clay-bearing  sediments,  where 
as  the  compressional  anisotropy  is  greater  than  the  shear  and  shear  orien 
tation  velocities  in  more  calcareous  sediments. 

Carlson  and  Christensen  (1979)  show  that  velocity  anisotropies 
are  well  developed  at  less  than  0.4  km  into  the  sediment  and  that  aniso¬ 
tropy  increases  markedly  with  depth  (depth  below  sea  level  is  not  meant 
here,  but  depth  of  overburden).  Concurrent  with  this  depth  dependent 
increase  in  anisotropies  are  related  increases  in  anisotropy  with 
increasing  density,  and  with  increasing  mean  velocity. 

Several  different  approaches  have  been  taken  toward  finding  a 
physical  mechanism  which  would  explain  the  existence  of  transversely  iso¬ 
tropic  sediments.  Hamilton  (1970)  considered  overburden  pressure  acting 
upon  clay  as  a  cause  of  such  a  system.  Since  clay  has  a  flocculated 
(collapsed  house  of  cards)  structure,  sufficient  pressure  should  lead  to 
an  ordering  of  the  clay  structure  in  which  the  platelets  would  assume 
positions  parallel  to  the  plane  of  bedding.  Hamilton  realized,  however, 
that  such  a  purely  mechanical  action  upon  clay  could  not  produce  the 
observed  anisotropies  at  depths  shallower  than  400-600  m,  although  a 
measurable  anisotropy  would  be  produced  at  these  depths.  This  explana¬ 
tion  of  sediment  velocity  anisotropy  has  not  proved  very  useful,  however, 
as  it  is  geared  to  explain  pure  silt  and  clay  sediments,  which  actually 
show  much  smaller  degrees  of  anisotropy  than  other  sediments;  a  non- 
calcareous  clay  from  0S0P  Site  210  shows  a  compressional  anisotropy  of 
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3.3%  and  a  shear  anisotropy  essentially  zero,  while  a  clay  enriched 
carbonate  ooze  and  chalk  from  the  same  site  shows  a  compressional 
anisotrophy  of  6.2%  and  a  shear  anisotrophy  of  6.63%. 

Carlson  and  Christensen  (1979)  proposed  mechanisms  for 
explaining  the  larger  velocity  anisotropies  in  calcareous  sediments.  They 
discounted  the  first  mechanism,  the  alignment  of  cracks  and  pores  in  the 
bedding  plane,  because  it  is  known  that  pores  and  cracks  close  under 
increasing  pressure;  this  would  cause  the  vertical  velocity  to  increase 
more  rapidly  than  horizontal  velocities  with  depth,  an  equality  (isotropy) 
eventually  being  reached.  Since  anisotropy  increases  with  depth,  this 
mechanism  cannot  be  the  driving  mechanism  for  transverse  isotropy. 

The  second  mechanism  examined  by  Carlson  and  Christensen  (1979) 
seem  to  be  the  most  adequate  in  explaining  the  observed  phenomena.  If 
the  cal  cite  crystals  in  calcareous  sediments  have  a  preferred  orientation, 
anisotropies  of  the  preferred  type  (i.e.,  transverse  isotropy)  would 
result.  Dandekar  (1968)  determined  that  cal  cite  crystals  exhibit  a 
velocity  of  5.61  km/ sec  parallel  to  their  c-axes,  and  a  velocity  of 
7.35  km/ sec  normal  to  their  c-axes  (the  c-axis  is  that  direction  along 
which  a  ray  of  light  is  not  doubly  refracted).  Therefore,  if  it  could  be 
ascertained  that  cal  cite  grains  align  their  c-axes  normal  to  the  plane  of 
bedding,  the  observed  anisotropies  would  be  accounted  for.  Manghnani, 

Schl anger,  and  Milholland  (1980)  reported  that  preliminary  X-ray  petro- 
fabric  data  for  carbonate  sediment  core  samples  from  DSDP  Site  289  show 
such  a  preferential  orientation  of  the  cal  cite  c-axes. 

Carlson  and  Christensen  (1979)  and  Manghnani,  Schlanger,  and 
Milholland  (1980)  pointed  out  that  calcareous  deep  sea  sediments  consist 
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almost  exclusively  of  biogenic  detritus.  Prominent  among  the  various 
detritus  in  calcareous  sediments  are  various  members  of  the  genus  Cocco- 
1 ithus  and  the  genus  Cyclococcolithina,  and  discoasters.  In  these  orga¬ 
nisms,  calcite  crystals  contained  in  the  disc  shaped  skeletons  are  aligned 
with  their  c-axes  normal  to  the  long  dimensions  of  the  microfossils. 

Though  these  disc-like  forms  are  randomly  oriented  in  oozes,  increasing 
compaction  leads  to  their  alignment  parallel  with  the  plane  of  bedding, 
leaving  the  calcite  crystal  c-axes  normal  to  the  plane  of  bedding. 

This  pressure  induced  alignment  of  c-axes  is  only  the  beginning 
of  the  ordering  of  carbonate  sediments,  however.  The  diagenetic  model  of 
Schlanger  and  Douglas  (1974)  posits  that  the  major  change  in  the  trans¬ 
formation  of  carbonate  ooze  to  chalk  to  limestone  is  the  solution  and 
reprecipitation  of  the  carbonate  as  cement.  Elimination  of  crystal 
defects  such  as  the  presence  of  strontium,  dissolution  of  small  crystals, 
and  the  growth  of  larger  crystals  reduce  the  free  energy  of  the  system. 
Furthermore,  Carlson  and  Christensen  (1979),  reporting  from  Bain  (1940) 
and  Bukry  (1971),  conclude  that  calcite  crystals  with  faces  nearly  per¬ 
pendicular  to  the  c-axes  are  more  resistant  to  dissolution  than  those  with 
faces  nearly  parallel.  Hence,  those  microforms  most  resistant  to  solution 
are  those  which  have  their  c-axes  normal  to  their  long  dimension.  In 
other  words,  the  process  of  solution  and  recrystallization  in  calcite 
operates  in  such  a  manner  that  most  calcite  c-axes  are  normal  to  the 
plane  of  bedding,  thus  producing  a  larger  velocity  of  propagation  parallel 
to  the  bedding  plane  than  normal  to  it. 

Milholland,  Manghnani ,  Schlanger,  and  Sutton  (1980)  believe 
that  the  role  of  clay  enrichment  in  increasing  velocity  anisotropy  is 
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best  explained  by  the  mechanism  proposed  in  Weyl  (1952).  According  to 
this  model,  clay  particles,  being  much  smaller  and  more  plastic  than  the 
cal  cite  crystals,  coat  cal  cite  grains  with  an  ooze,  reducing  the  point- 
to-point  contacts  of  crystals,  and  increasing  the  total  reactive  surface, 
thus  catalyzing  the  solution  and  recrystallization  of  cal  cite  at  a  rate 
more  rapid  than  normal. 

In  sunmary,  many  measurements  of  anisotropy  establish  values  of 
0-15%  with  10%  being  typical.  Many  mechanisms  have  been  proposed  to 
account  for  such  anisotropies,  including  pore  and  crack  alignment,  over¬ 
burden  pressure,  and  the  alignment  of  cal  cite  crystals.  We  conclude 
that,  while  overburden  pressure  plays  a  minor  role  in  determining  veloc¬ 
ity  anisotropies,  the  principal  cause  of  such  anisotropies  is  the  pre¬ 
sence  of  calcite  crystals  aligned  with  their  c-axes  normal  to  the  plane 
of  bedding.  The  mechanism  that  best  explains  this  is  the  deposition  of 
certain  biogenic  detritus  and  their  subsequent  dissolution  and  re- 
crystallization  with  the  proper  alignment.  The  existence  of  such  align¬ 
ments  has  been  confirmed  experimentally. 
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CHAPTER  III 


THEORY  OF  REFLECTION  FROM  TRANSVERSELY  ANISOTROPIC  SOLIDS 

From  the  data  and  geological  models  available,  it  is  reasonable 
to  model  an  anisotropic  sediment  as  a  transversely  isotropic  solid,  that 
is,  as  a  body  containing  a  single  unique  axis,  all  rotations  about  the 
axis  being  equivalent.  This  axis  will  be  taken  to  be  the  z-axis,  per¬ 
pendicular  to  the  plane  of  bedding. 

A.  Independent  Elastic  Parameters  for  Transverse  Isotropy 

The  first  problem  which  arises  from  a  departure  from  isotropy 
is  an  increase  in  the  number  of  independent  elastic  parameters  of  the 
system.  While  an  isotropic  system  has  but  two  independent  elastic  param¬ 
eters,  a  transversely  anisotropic  system  has  five.  In  obtaining  these 
constants,  we  shall  follow  the  development  of  Love  (1944),  in  a  somewhat 
more  modern  notation.  If  one  begins  from  the  stress-strain  relation, 
Hooke's  Law, 


Sij  =  Cijklekl 

there  are  a  total  of  81  possible  elastic  parameters.  Here  S^.,  e^ ,  and 
Sjkl  are  comPonents  the  stress,  strain,  and  elastic  tensors, 
respectively,  and  repeated  indices  are  summed  over.  However,  a  set  of 
equivalences. 


ppy.BTYiMft.  pjtOS  BiHHKH  HU® 


Ci jkl  =  Ckli j  =  Cjikl  =  Cijlk 

reduces  the  number  of  independent  elastic  parameters  to  21  in  the  case  of 
complete  anisotropy.  The  set  of  equivalences  springs,  respectively,  from 
the  conditions  that  a  thermodynamic  internal  energy  exists,  that  the 
stress  S.jj  is  symmetric,  and  that  the  strain  e^  is  symmetric  under  inter¬ 
change  of  index. 

To  simplify  the  notation,  the  elements  of  these  tensors  are 
identified  with  the  components  of  a  symmetric  6x6  matrix.  The  fol¬ 
lowing  associations  of  subscripts  are  made:  xx=l ,  yy=2,  zz=3,  yz=4, 
zx=5,  xy=6.  Thus,  the  stress-strain  relation  may  be  rewritten  using 
this  notation  as 


SI  =  CIJej 

To  further  reduce  the  number  of  independent  elastic  parameters, 
let  us  utilize  the  strain-energy  function,  defined  by  the  relation 

c  _  JW_ 

*1  "  aej 


As  discussed  in  Love  (1944)  and  Aki  and  Richards  (1980)  the  actual  quan¬ 
tity  chosen  for  W  may  be  either  the  free  energy  or  the  internal  energy, 
depending  upon  whether  conditions  are  isothermal  or  adiabatic.  It  is 
standard  in  most  texts  to  choose  adiabatic  conditions,  since  at 
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wavelengths  greater  than  a  few  millimeters  the  thermal  diffusion 
coefficient  in  sediments  is  much  longer  than  the  period  of  elastic  waves. 
Combining  Hooke's  Law  with  the  defining  relation  for  W,  one  obtains  the 
result: 


2W  =  CIJeIeJ 


Written  in  its  entirety,  this  equation  is: 


2W  =  C1-je1  +  2C^^2  +  2C13elC3  +  2C14ele4 


+  ZC15ele5  +  2C1 6ele6  +  C22e2  +  2C23e2e3 


+  2C24e2e4  +  2C25e2e5  +  2C2ge2eg  +  C33e3 


+  2C34e3e4  +  2C35e3eg  +  2C36e3eg  +  C44e4 


+  2C45e4e5  +  2C46e4e6  +  C55e5 


+  2C56e5e6  +  C66e6 


We  now  demand,  for  transverse  anisotropy,  that  this  quantity  be  invariant 
under  a  rotation  about  the  z-axis.  Under  a  rotation  of  angle  0  about  the 
z-axis,  the  strains  take  the  form 
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2  2 

g 1 i  =  e^cos  e  +  e2sin  e  +  egsmecose  , 

.2  2 

e2,  =  e^sin  e  +  e2cos  e  -  egSinecose  , 
e3’  =  e3  » 

e^,  =  e^cose  -  e5sine  , 

e5,  =  e^sine  +  e5cose  , 

2  2 

e6,  =  —  2e i si necose  +  2e2sinecose  +  e6(cos  e-sin  e)  , 

where  primes  refer  to  the  values  of  the  strains  in  the  rotated  system. 

If  these  substitutions  are  made  into  the  strain-energy  function 
and  coefficients  of  corresponding  strain  products  are  equated,  a  set  of 
21  equations  is  obtained.  Fortunately,  this  set  of  equations  can  be  par¬ 
titioned  into  smaller  subsets,  corresponding  to  a  partition  of  the  set  of 
elastic  parameters. 

Consider  the  subset  obtained  by  equating  the  coefficients  of 
the  products  e^,  e2e2,  e^e2,  e^g,  e2eg,  and  egeg,  respectively: 


4  2  2  3 

C-j ^  =  C-|-]COS  e  +  2C12cos  esin  e  -  4C-jgSin9cos  e 


4  3  2 

+  C22sin  e  -  4C2gSin  9cos9  +  4CggSin  9cose  , 
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4  2  2 

C22  =  cnsin  6  +  2C12cos  esin  0  +  4C-|gSin  ecose 


4  3  2  2 

+  C22COS  0  +  4C2gSin0cos  0  +  4Cgg sin  ©cos  0  , 


2  2  4  4  2  2 
C] 2  =  C-^cos  0sin  0  +  C1?(cos  0+sin  e)  +  2Clfisin0cose(cos  0-sin  0) 


'12' 


'16- 


2  2  2  2  2  2 
+  C22Sin  0cos  0  -  2C2gSin0cos0(cos  0-sin  0)  -  4CggSin  0cos  0  , 


3  2  2  4  2  2 

C16  =  C^sinocos  0  +  C^sinocosoCsin  0-cos  0)  +  C^g(cos  0-3cos  0sin  0) 


3  4  2  2  2  2 

-  C22sin  0cos0  -  c2g (sin  0-3cos  0sin  0)  -  2CggSin0cos0(cos  0-sin  0)  , 


'26' 


3  2  2  4  2  2 

C26  =  cnsin  0cose  +  C^2Sin0cos0(cos  0-sin  0)  -  C^g(sin  0-3cos  0sin  0) 


3  4  2  2  2  2 

-  C22sin0cos  0  +  C9A(cos  0-3cos  osin  0)  +  2Cfifisin0cos0(cos  0-sin  0)  , 


'26' 


'66' 


2  2  2  2  2  2 
C66  =  cnsin  0COS  9  "  ZC^sin  0cos  0  +  4C^gSin0cos0 (cos  0-sin  0) 


2  2  2  2  2  2  2 
+  C22Sin  0cos  0  -  4C2gSin0cos0(cos  0-sin  0)  +  Cgg(cos  0-sin  0) 


Adding  the  fifth  and  sixth  equations,  then  subtracting  the  second  equation 
from  the  first,  and  removing  a  factor  of  sin©  from  both  results  yields 

0  =  (C^ ^ -C22 ) C°S6  -  2(C1g+C2g)sin0  , 

0  =  (cy|-C22^cose  +  2(C^g+C2g)sin0 
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Solving  for  C^-C^  and  C^+C^  yields 

C11  =  C22 
C26  =  C16 


Substituting  these  results  into  the  third  and  sixth  equations  and  com¬ 
bining  like  terms,  one  obtains 


2  2  2  2 
0  =  2C^-jCos  esin  e  -  2C-j2sin  ecos  e 


2  2  2  2 
+  4C.|gCOs8sin0(cos  e-sin  e)  -  4Cggsin  ecos  e  , 


0  =  2C11cos2esin2e  -  2C12sin2ecos2e 

2  2  2  2 
+  8CigCosesine(cos  e-sin  e)  -  4Cggsin  ecos  e 


On  substituting  these  two  equations,  it  is  obvious  that 


2  2 

Making  this  substitution  and  removing  a  factor  of  2cos  esin  e,  one 
obtains 


C66  =  V2(CirC12) 
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Let  us  next  consider  the  second  set  of  equations,  obtained  by 
equating  the  coefficients  by  e^,  e2e3»  an<*  e3e6*  respectively > 


2  2 

C-J3  c  C13cos  0  +  ^23Sin  0  -  2C3gSinecose  , 


2  2 

C23  =  C^sin  0  +  C23COS  e  +  2C3gSinecos0  , 


2  2 

C36  =  C-|3sinecose  -  CgjSinecose  +  C3g(cos  e-sin  0) 


Simplifying  the  first  two  equations  yields 


0  =  (c23-c  ^sinS  -  2C3gSin8cose  , 


0  =  (c13-c  23)cos^0  -  2C3gSin0cose  , 


which  are  correct  at  sine=0  or  cose=0  only  if 


C13  "  C23 


from  which  it  follows 


C36  -  0 


The  third  set  of  equations,  arising  from  equating  coefficients 


of  ,  e^€g ,  and  egeg ,  is 
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2  2 

c44  =  C44cos  6  +  2C45cosesine  +  C55sin  e  , 

2  2 

C45  =  C^cosesine  +  C^g(cos  e-sin  e)  +  Cggsinecose  , 

2  2 

C55  =  C44s1n  0  “  2C45sinecose  +  CggCos  e 


Simplifying  the  first  two  equations  yields 

0  =  -Cosine  +  2C45cos6  +  Cosine 
0  =  -C^cose  -  2C^gSine  +  C55cose 


which  has  a  solution  consistent  with  the  third  equation. 


C44  =  C55 


and 


The  fourth  set  of  equations,  obtained  from  equating  coefficients 
by  e3e4  and  e3e5,  is 

C34  =  C34cos9  +  C35sine  » 
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:*-.*•* 


C35  =  -C34sine  +  C35cose 


Consistency  at  cose=0  demands 

C34  =  C35  =  0 


The  fifth  subset  of  equations  yields  no  new  information.  The  final  sub¬ 
set,  obtained  by  equating  coefficients  by  e^,  e-jeg,  e^,  e2e5*  e4e6* 
and  e5e6,  yields: 

3  2  2 

C-J4  =  C^cos  e  +  C15sin0cos  e  +  c^sin  ecose 

3  2  3 

+  C25sin  e  -  2C46sin0cos  0  -  2CggSin  0cos0  , 

2  3  3 

C15  =  ■C'|4sin0cos  0  +  C15cos  0  -  C24sin  0 

2  2  2 
+  C25sin  0cos0  +  2C4gSin  0cos0  -  2CggSinecos  0  , 

2  3  3 

C24  =  ci4sin  0COS0  +  ci5sin  9  +  C^cos  0 

2  .2  2 
+  C2gSin0cos  0  +  2C^gSin0cos  0  +  2CggSin  0cos0  , 

.32  2 

^25  =  "('l4sin  9  +  ^I5sin  9cose  ~  C24sin0cos  0 

3  2  2 

+  C25cos  0  -  2C4gSin  0COS0  +  2C5gsin0cos  0  , 
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which  can  be  consistent  only  if 


'1 5 


C24  "  0 


Combining  all  these  results,  the  strain-energy  function  may, 
for  transversely  isotropic  media,  be  written  in  terms  of  only  five  param¬ 
eters  (^ii »  ^33 »  ^13*  C44) * 


2W  -  C11(e1+e2)  +  C33e3  +  2C-|  3(e]+e2)e3 


+  2Ci 2eie2  +  c44(e4+e5)  +  c66e6 


where  Cgg  =  1  /2 ( C -j  1  -C-^) ♦ 

B.  Equations  of  Motion  for  Transversely  Isotropic  Media 

Reverting  to  our  initial  notation,  the  equation  of  motion  may 
be  obtained  from  the  potential  energy  via  the  stresses: 


where  S- •  =  3w/3e, . ,  p  is  the  density,  the  quantities  u-  are  the  com- 

I  J  I  J 

ponents  of  displacement,  and  t  is  the  time.  Substituting 


and 
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3U  •  3U  . 

e .  .  -  — +  — 4- 

1J  3XJ  3X1 


into  the  equations  of  motion,  one  obtains  the  set  of  equations 


32u 


3t 


_3_ 

3X 


r  3U.  +  r  3V 
3x  +  i 2  3y 


+ 


3W 
13  3Z 


where  u,  v,  w  have  been  substituted  for  the  components  of  the  displace¬ 
ment  vectors  u  ,  u  ,  u  ,  respectively. 
a  y  z 

In  order  to  attach  a  bit  more  physical  intuition  to  these 
equations,  let  us  examine  waves  whose  propagation  and  polarization  are 
both  along  the  coordinate  axes. 
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For  an  x-propagating  (v=w=0)  or  a  y-propagating  (u=w=0) 
compressional  wave,  the  equations  of  motion  reduce  to 


p 


or 


32V  _  c 


a2v 

n  ^7 


This  identifies  C^/p  as  the  square  of  the  horizontal  compressional 
velocity. 

For  a  z-propagating  (u*v=0)  compressional  wave,  the  equations 

become 


2  2 
3  W  _  r  3  U 

P  ,  C33  2 


3t 


3Z 


This  identifies  C33/P  as  the  square  of  the  vertical  compressional 
velocity. 

For  shear  waves  which  either  propagate  in,  or  are  polarized  in 
the  z  direction,  the  differential  equations  reduce  to 
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This  identifies  C^/p  as  the  square  of  the  shear  velocity  in  the  vertical 
plane. 

Shear  waves  which  propagate  in,  and  are  polarized  in,  the  plane 
of  bedding  obey 


2  2 
3  U  _  r  3  U 

'  C66  77 

3t  3y 


* 


32V  p  3^V 

P  at2  ‘  66  ax2 


This  identifies  Cgg/p  as  the  square  of  the  shear  velocity  in  the  plane 
of  bedding. 

In  summary,  four  of  the  five  independent  elastic  parameters  of 
an  anisotropic  medium  may  be  associated  with  compressional  and  shear 
velocities  in  axial  directions,  as  follows: 
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In  the  vertical  plane,  the  constant  Cgg  does  not  play  a  role, 
as  it  appears  only  for  motion  entirely  in  the  plane  of  bedding.  There¬ 
fore,  our  first  important  result  is  that  shear  anisotropy  does  not  play 
a  role  in  the  problem.  The  anisotropic  properties  of  our  tranversely 
anisotropic  medium  affect  motion  only  through  the  compressional  anisotropy 
and  in  the  behavior  of  the  constant  C^. 

Let  us  assume  a  plane  wave  form  for  the  displacements, 

^  exp(iqz+ikxx-iwt)  , 

where  kx  is  the  horizontal  wave  number,  assumed  from  measurements 
in  water,  q  is  the  vertical  wave  number,  still  to  be  identified,  and  u> 
is  the  angular  frequency. 

Substitution  of  this  form  into  the  equations  of  motion  yields 

2  2  2 

-pa)  u0  =  -Cnkxu0  -  C44q  u0  -  (C13+C44)  kxqw0 

-P“2w0  =  'C44kxwo  '  C33q2wo  "  ^C13+C44^  kxquo  • 

To  simplify  these  equations,  define 


These  substitutions  yield 


S  = 


2  r  k2 
pu>  -  C44kx 


^T-C44q  ^  uo  "  Jqkxwo  =  0  ’ 

(S‘C33q2)  wo  "  Jqkxuo  =  0  ’ 

which  have  solution  only  if 

(r-C44q2)(s-C33q2^  -  Ay  =  0  . 

This  consistency  equation  defines  the  allowed  vertical  wave  numbers.  It 
2 

is  a  quadratic  in  q  ,  with  solutions 


r  + 


/r2  ■  4C33C44TS 
2C33C44 


and 


4C33C44TS 


2C33C44 
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where 


r  =  TC33  +  SC44  +  J  kx 

Thus,  the  plane  wave  solution  to  the  equations  of  motion  has 
two  vertical  waves  allowed  and  takes  the  form 


exp(iq1z) 


+ 


exp(iq2z). 


exp(ikxx-iwt) 


In  order  to  check  this  solution  for  the  vertical  wave  number, 
let  us  examine  the  isotropic  limit  by  setting  C^=C33=x+2u,  C13=X,  and 
C44=VI»  where  x  and  y  are  the  usual  Lame  parameters  for  isotropic  media. 
The  wave  numbers  for  the  two  allowed  plane  waves  obtained  by  such  a  sub¬ 
stitution  are: 


2  2  2 
q-j  =  (pa  -ykx)/y  , 

q 2  =  (pu>^-(x+2y ) )/y 

Thus  the  plane  wave,  in  the  isotropic  case,  becomes  the  familiar  purely 

1  /2 

congressional  wave  of  velocity  ((x+2u)/p)  '  ,  and  purely  shear  wave  of 
1  /2 

velocity  (y/p)  .  Thus,  q-j  is  associated  with  a  quasi-shear  wave,  whil 

q2  is  associated  with  a  quasi -compressional  wave. 


28 


This  development  has  already  revealed  several  points  of 
interest.  The  shear  anisotropy  does  not  enter  the  problem  of  the  propaga¬ 
tion  of  elastic  waves  in  the  vertical  plane.  The  coupled  wave  equations 
for  the  displacements  differ  from  the  isotropic  equations  only  by  having 
C-j-j  and  C33  where  the  isotropic  equations  would  have  a  single  constant. 
This  difference  is  sufficient  to  prevent  a  useful  separation  in  poten¬ 
tials;  that  is,  purely  shear  and  compressional  waves  cannot  be  defined. 

D.  Reflection  of  a  Plane  Wave  from  an  Anisotropic  Solid 

Now  that  the  behavior  of  plane  waves  in  a  single  infinite 
medium  has  been  examined,  let  us  proceed  to  derive  the  reflection  co¬ 
efficient  for  plane  waves  from  the  fluid  side  incident  upon  the  boundary 
between  a  semi -infinite  isotropic  ocean  and  a  semi -infinite,  transversely 
anisotropic,  solid  ocean  floor.  The  coordinate  system  in  which  we  shall 
work  is  standard,  the  x-y  plane  being  the  boundary  plane  and  the  positive 
z-axis  extending  down  into  the  solid. 

The  standard  boundary  conditions  apply;  that  is,  vertical 
displacement,  vertical  stress,  and  tangential  stress  are  continuous 
across  the  interface.  Since  the  liquid  layer  cannot  support  shear  waves, 
the  tangential  stress  must  vanish  at  the  interface.  The  boundary  condi¬ 
tions  may  then  be  written  in  the  form: 

- ( 1 -R) i kQ  =  w1  +  w2 

(l+R)PgW  =  i  C22  (d'j  +q2w2 )  +  iC-)  3^  ^2  ^  ’ 

0  =  q-jU-j  +  q2u2  +  kx(w1+w2) 
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where  R  is  the  reflection  coefficient,  pQ  is  the  water  density,  and  kQ  is 
the  vertical  wave  number  in  the  water.  These  equations  correspond, 
respectively,  to  continuity  of  vertical  displacement,  vertical  stress, 
and  horizontal  stress.  There  are  five  unknowns  (R,  u^ ,  U2>  w^ ,  w^),  and 
only  three  boundary  equations.  The  equations  of  motion  provide  the 
remaining  two  equations  needed  to  solve  for  the  unknowns.  Through  the 
relations  obtained  before  we  have: 

0  *  (T-C44ql)ul  '  JkxVl  ■ 

0  =  (S‘C33q!)w2  '  Jkxq2u2 


The  system  of  equations  is  now  soluble,  and  the  reflection  coefficient 
is  given  by  the  ratio  of  determinants: 
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'ipoai  C33ql  C33q2 


0  k. 


0  -Jkxql 


C]3kx 


T-C44"? 


C13kx 


S-C33q2 


_Jkxq 


2 

ipo“  C33ql  C33q2  C13kx  C13kx 


0  k. 


0  -Jk  q,  0  T-C..q., 


0  ^23^2 


-Jkxq2 


Upon  working  the  algebra  entailed  by  these  determinants,  one  obtains 


R  -  ft  -  B 
K  7T+T 


where 


A  =  k 


o  Jkx^q2_ql^  C33q2 (T-C44ql^  '  C1 3ql  (S_C33ql) 


(C33qlq2‘C13kx)Jj2|<xqlq2  ‘  (S-C33ql)(T_C44q7  )| 


;;  v«**c 
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and 


-q2(T'C44ql)(Jkx+S-C33q2)  +  Jkxq2(T+C13qlq2 


On  examining  the  isotropic  media  limit  of  this  result  (making 
the  same  substitutions  made  in  the  reduction  of  the  wave  number),  one 
obtains 


which  is  the  usual  form  for  the  reflection  coefficient  found  in  Ewing, 
Jardetsky,  and  Press  (1957). 

Thus,  we  have  found  that  we  can  solve  for  the  reflection 
coefficient,  applying  standard  matrix  techniques  to  a  system  of  equations 
containing  the  standard  boundary  conditions  and  two  additional  equations 
which  our  initial  demand  upon  the  form  of  the  solution  forced  from  the 
differential  equations.  Our  solution  is  somewhat  complicated,  but  reduces 
in  the  isotropic  limit  to  a  familiar  form. 

E.  Reflection  from  a  Layered  Anisotropic  Solid 

The  reflection  coefficient  having  been  obtained  for  two 
semi-infinite  media,  the  next  step  is  to  expand  the  theoretical  frame¬ 
work  to  include  finite  layers. 
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In  a  finite  layer,  the  displacements  now  include  upward  and 


downward  propagating  plane  waves  and  take  the  form 

w  =  w-j  exp(iq^z)  +  w2  exp(-iq-jz) 

+  w3  exp(iq2z)  +  w4  exp(-iq2z) 

u  =  u1  exp(iq-| z)  +  u2  exp(-iq2z) 

+  u3  exp(iq2z)  +  u4  exp(-iq2z) 

from  which  we  may  obtain  the  normal  stress 

Szz  =  ^C13kxul+C33qlwl^  exP^ 1  2 ) 

^  ^ '^1 3kx^2_^33ql ^2 ^ 

+  1 (Ci3kxu3+C33q2w3)  exp(iq2z) 

+  1<c13kxu4-C33q2w4)  “Pl-’V1 
and  the  tangential  stress 

Szx  =  1  WVWl)  exp(iqlz> 

+  i  c44(kxw2“qiu2^  expC-iq] z) 


+  i  c44(kxw3+q2u3^  exP^1c^22^ 

+  i  C44(kxw4-q2u4)  exp(-iq2z) 

We  also  have  available  the  relations  between  the  vertical  and 
horizontal  displacements  which  must  be  satisfied  for  the  existence  of  a 
plane  wave  solution.  These  relations  are  similar  to  those  used  earlier 
to  solve  the  problem  of  two  semi-infinite  media.  They  are: 

0  =  (T-C44ql)Ul  "  Jqlkxwl 

0  =  (T-C44qi )u2  +  Jqlkxw2  ’ 

0  =  (S_C33q2)u3  "  Jq2kxu3  * 

0  =  (S_C33q2)u4  +  Jq2kxu4 

These  consistency  relations  can  be  used  to  eliminate  the 
horizontal  displacements  (u^ )  in  terms  of  the  vertical  displacement  in 
the  stresses  Szz  and  Sxz-  This  done,  we  define  a  vector  of  quantities 
which  are  continuous  across  an  interface  between  two  transversely 
anisotropic  solids. 

Y  = 


34 


This  vector  can  now  be  written  entirely  in  terms  of  the  vertical 
displacements  (w^ )  by  using  a  matrix  format. 


Y  = 


x 


1111 
^C13kxXl+C33ql  ^  '  1^C13KxXl+C33ql^  ^C13kxX2+C33q2^  "  1  ^C13kxX2+C33q2^ 

lC44^kx+Xl ql )  lC44^kx+Xlql ^  1 C44^kx+X2q2^  lC44^kx+X2q2^ 

X1  -X1  X2  "  X2 


"g 

0 

0 

0 

W1 

0 

1/G 

0 

0 

W2 

0 

0 

Q 

0 

w3 

0 

0 

0 

1/Q 

*4. 

or  Y=E  D(z)B,  where  E  is  the  first  matrix,  D(z)  is  the  second  matrix,  and 
B  is  the  vector  containing  the  vertical  displacements.  The  quantities  in 
the  matrices  are: 

G  =  exp(iq-|Z) 

Q  =  exp(iq2z) 

X1  = 

X2  =  (S_C33qi)/Jq2kx  * 

The  vector  Y  has  two  interesting  properties.  First,  it 
embodies  the  solutions  to  the  differential  equations  for  plane  wave 
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motion  in  a  transversely  anisotropic  media.  It  thus  expresses  the  values 
of  displacement  and  stress  at  each  point  z  within  a  layer.  Sfecond,  it 
contains  only  quantities  that  are  continuous  across  the  interface  between 
two  solids. 

Using  these  properties  we  can  now  develop  the  method  for 

computing  the  reflection  coefficient  of  a  plane  wave  in  the  ocean  incident 

upon  a  homogeneously  layered,  transversely  isotropic  ocean  floor.  Let 

each  layer  be  numbered  ^rom  1  to  n  with  the  ocean  being  layer  0  and  the 

semi-infinite  substrate  being  layer  n+1 .  At  the  top  of  the  layer  labeled 

m,  the  depth  dependent  quantities  are  represented  as  Ym(0),  and  at  the 

bottom  of  the  layer,  as  Y  (Hi,  where  H  is  the  thickness  of  the  layer. 

m  m  m 

Thus,  the  continuity  of  the  stress  and  displacement  across  the  inter¬ 
face  between  layers  m  and  m+1  is  given  by  Y  ^ (0)=Ym(H)  or 

Em*l  Vl  ‘  Wm  • 


This  can  be  rewritten  as 


®m+l  •  Cl  W'm  • 

where  E~j-|  is  the  inverse  of  the  matrix  Em+i . 

In  the  substrate  layer,  the  fact  that  no  further  reflections 
occur  reduces  the  displacements  to  downward  propagating  plane  waves  only: 

W  =  A-j  exp(iq^z)  +  A^  exp(iq2z) 

U  =  A3  exp{iq1z)  +  A4  exp(iq2z) 


Ifc. 
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The  consistency  relations  which  allow  plane  wave  solutions  are 


(T-C44qq)Al  '  JkxqlA3 

(s-c33q2)A2  -  Jkxq2A4 

yielding,  upon  substitution  into  the  definition  of  stress,  the  vector 


’1 

rl 

<C33ql+C13kxV 

(C33q2+C13kxX2^ 

-< 

+ 

II 

iC44^kx+qlXl ^ 

+  A2 

iC44(kx+q2X2) 

*1 

h 

where  N  and  L  are  the  first  and  second  vectors,  respectively,  in  the 
equation  for  above,  and 


X1  =  (T'c44q?)/Jkxql 
X2  =  (S'C33q2)/Jkxq2 


Therefore,  in  terms  of  the  amplitudes  at  the  top  of  the  first  solid 
layer,  the  substrate  stresses  and  displacements  are 


Yn+i  =  A]N  +  A2L  - 


1 


TT ( E  D  ( H  ) E” 1 ) 
„  ,  m  m'  m  m 
m=l 


Yn(0) 


=  PB, 


where  P  is  the  product  matrix  through  which  the  unknown  amplitudes  in  the 
intermediate  layers  are  eliminated. 

In  the  top  semi-infinite  fluid  layer,  however,  we  have 


where 


‘-ik0(l-R) 

p0w2(HR) 

0 

0 


and 
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Y1  - 
o 


ik 


pow 


Since  the  fluid  does  not  support  tangential  displacements,  the  transition 
of  the  first  interface  is  handled  by  writing 


v°>*v 


where  uQ  is  the  unconstrained  horizontal  displacement  at  the  top  of  the 
first  layer.  Substituting  for  Y^(0)  in  the  equation  for  Yn+1  yields  the 
system  of  linear  equations  which  must  be  solved  to  obtain  the  reflection 
coefficient 


V  +  A2L  -  PViR  -  P-4Uo  =  PVo  • 

where  P.^  is  the  fourth  column  of  the  matrix  P.  As  there  are  only  four 
unknown  quantities  in  this  relation  (A^  ,  A^,  R,  uQ),  a  solution  may  be 
obtained  via  any  number  of  standard  methods  for  linear  equations. 

Thus,  we  have  established  a  method  by  which  the  reflection 
coefficient  can  be  obtained  for  a  system  containing  an  arbitrary  number 


39 


of  layers  in  which  the  unknown  displacement  amplitude  of  intermediate 
layers  do  not  appear.  This  method  differs  from  that  for  isotropic 
layers  in  that  equations  relating  to  displacement  amplitude  obtained  from 
the  differential  equations  must  be  used  to  reduce  the  matrices  to  a  form 
in  which  the  movement  from  one  layer  to  another  can  be  handled  by  a  sim¬ 
ple  multiplication  of  matrices.  The  evaluation  of  the  formal  solution  is 
obviously  a  matter  for  the  computer. 
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CHAPTER  IV 


PROPERTIES  OF  THE  REFLECTION  COEFFICIENTS 

In  our  computation  of  reflection  coefficients  we  have  utilized 
a  standard  environment.  The  semi -infinite  layer  in  which  the  plane 
waves  originate  is  water,  with  a  sound  speed  of  1550  m/sec  and  a  density 
of  1.043  g/cm  .  The  underlying  solid  layers  have  elastic  parameters 
which  are  discussed  below.  Our  computations  have  been  made  at  a  standard 
frequency  of  50  Hz. 

The  primary  difficulty  in  obtaining  meaningful  results  from 
the  theory  developed  above  lies  in  the  paucity  of  published  data  on 
values  for  the  elastic  parameters.  We  have  located,  to  date,  only  three 
full  sets  of  elastic  parameters.  The  earliest  were  obtained  from  Podio 
et  al.  (1968),  and  are  recapitulated  in  Tables  II  and  III.  These  two 
samples  correspond,  respectively,  to  compressional  anisotropies  (computed 
according  to  the  formula  in  Chapter  II)  of  2.81%  and  3.89%.  The  latest 
set  of  parameters  were  obtained  by  Jones  and  Wang  (1981)  for  cretaceous 
shales  from  the  Williston  Basin,  and  are  recapitulated  in  Table  IV. 

This  set  of  parameters  corresponds  to  a  compressional  anisotropy  of 
20.6%.  Thus,  while  our  samples  are  few,  nevertheless  realms  of  both 
high  and  of  low  anisotropy  are  represented. 

In  Section  IV. A,  we  discuss  the  effects  sediment  anisotropy 
has  on  the  vertical  wave  number,  and  note  the  possibility  of  backward 
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TABLE  II 


ANISOTROPIC  ELASTIC  PARAMETERS  FOR  GREEN  RIVER  SHALE, 
MEASURED  UNDER  8000  psi  CONFINING  PRESSURE* 

Shale  Density  2.35  g/cm^ 


Parameter 

Value 

(kbars) 

Velocity 

(m/sec) 

C11 

476 

4500 

C33 

450 

4243 

C44 

154 

2556 

C13 

170 

2690 

♦Data  from  Podio  et  al .  (1968) 
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TABLE  III 

ANISOTROPIC  ELASTIC  PARAMETERS  FOR  GREEN  RIVER  SHALE 
MEASURED  UNDER  1000  psi  CONFINING  PRESSURE* 
Shale  Density  2.35  g/cm^ 


Parameter 


Value 

kbars 


Velocity 

(m/sec) 


TABLE  IV 


ANISOTROPIC  ELASTIC  PARAMETERS  FOR 
WILLISTON  BASIN  CRETACEOUS  SHALE* 
Shale  Density  2.42  g/cm 


Parameter 

Value 

(kbars) 

Velocity 

(m/sec) 

C11 

343 

3765 

C33 

227 

3063 

C44 

54 

1494 

C13 

107 

2103 

*Data  from  Jones  and  Wang  (1981) 
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traveling  phase  fronts.  Section  IV. B  deals  with  the  effects  of 
anisotropy  on  reflection  from  hard  sediments.  In  Section  IV. C  a  scheme  to 
approximate  the  constant  is  considered,  and  in  Section  IV. D,  the 
effects  of  anisotropy  on  reflection  from  soft  sediments  is  examined. 

A.  Effects  of  Anisotropy  on  Wave  Numbers 

Since  the  vertical  wave  numbers  and  q2  only  reduce  to  their 
normal  form  in  the  limit  of  isotropy,  it  is  of  interest  to  observe 
their  behavior  for  the  transversely  anisotropic  media  before  proceeding 
to  the  computation  of  the  surface  reflection  coefficient. 

Plots  in  the  complex  plane  of  q^  and  q2  for  Green  River  Shale 
under  8000  psi  confining  pressure  are  contained  in  Fig.  1.  In  this 
figure,  solid  lines  represent  regions  of  overlays  of  q^  and  q2,  dotted 
lines  represent  q1  only,  and  dashed  lines  represent  q2  only;  the  grazing 
angles  are  marked  for  reference.  Recall  that  q1  is  the  quasi -shear  ver- 
tical  wave  number  and  q  is  the  quasi -compressional  vertical  wave  number, 
based  on  their  behavior  in  the  isotropic  limit.  Figure  1  shows  that,  as 
in  the  isotropic  case,  each  wave  number  has  a  purely  real  part  cor¬ 
responding  to  transmission  into  the  sediment  and  a  purely  imaginary  part 
corresponding  to  total  reflection  from  the  sediment.  The  total  reflec¬ 
tion  occurs  due  to  velocities  in  sediment  being  greater  than  that  in 
water.  The  anisotropic  case  also  has  a  region  at  very  low  grazing  angle 
(in  this  case  up  to  10°)  in  which  the  vertical  wave  numbers  are  neither 
purely  real  nor  purely  imaginary.  This  region  does  not  occur  for  iso¬ 
tropic  media.  In  this  range  q2  is  approximately  equal  to  -q^. 

The  existence  of  negative  values  for  the  real  part  of  the 
vertical  wave  number  q-j  indicates  that  in  this  range  there  are  backward 
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FIGURE  1 

THE  VERTICAL  WAVE  NUMBERS  IN  THE  COMPLEX  PLANE 
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traveling  phase  fronts,  i.e.,  phase  fronts  traveling  toward  the  interface 
from  the  solid  side.  As  will  be  discussed  below,  this  quirk  of  the  ver¬ 
tical  wave  number  does  not  affect  the  reflection  coefficient. 

The  off-axis  complex  vertical  wave  numbers  occurred  for  input 
data  in  which  the  velonty  contrast  between  media  was  very  high.  For 
Green  River  shale  at  1000  psi  confining  pressure,  the  region  in  which  the 
vertical  wave  numbers  are  neither  purely  real  nor  purely  imaginary  becomes 
much  shorter,  and  for  Willi ston  Basin  shale,  which  has  no  zone  of  total 
reflection,  it  vanishes  altogether.  Thus,  such  a  phenomenon  would  most 
likely  occur  in  an  oceanic  environment  of  exposed  basalts  or  hard  sedi¬ 
ments  rather  than  in  the  normal  gradual  layering  from  ooze  to  hard  sedi¬ 
ments  to  basalt. 

Two  important  results  appear  from  this  study  of  the  vertical 
wave  number.  First,  in  hard  sediments  we  may  see,  at  low  grazing  angles, 
wave  numbers  which  represent  propagating  but  heavily  damped  waves.  Sec¬ 
ond,  waves  of  this  type,  which  corresponds  to  a  quasi -shear  wave,  will 
be  propagated  backwards,  towards  the  interface. 

B.  Effects  of  Anisotropy  on  Fluid-Solid  Interface  Reflection 
Coefficient 

In  computing  the  surface  reflection  coefficient  R,  we  will  find 
that  measured  anisotropies  produce  observable  changes  in  both  the  magni¬ 
tude  and  phase  of  R  when  compared  to  the  isotropic  case.  We  will  compute 
bottom  loss  curves  for  such  surfaces,  where  BL  =  -20  log-jQ  |R|.  The 
sediments  studied  here  are  those  for  which  full  sets  of  anisotropic 
parameters  have  been  measured,  i.e.,  the  shales  described  in  Tables  II- 
IV.  These  materials  do  not  correspond  to  the  soft,  unconsolidated 
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material  found  over  much  of  the  ocean  floor,  but  are  more  typical  of 
exposed  basalt  and  outcropping  sedimentary  rocks. 

In  Figs.  2-7,  we  compare  computations  of  bottom  loss  and  the 
phase  of  the  reflection  coefficient  for  isotropic  and  anisotropic  con¬ 
solidated  sediments.  On  each  figure  solid  curves  represent  isotropic 
computations ,  and  dashed  curves  represent  anisotropic  curves.  Figures  2 
and  3  show  bottom  loss  and  phase  of  the  reflection  coefficient,  respec¬ 
tively,  for  a  Green  River  shale  under  8000  psi  confining  pressure; 

Figs.  4  and  5  show  these  quantities  for  a  Green  River  shale  under 
1000  psi  confining  pressure;  and  Figs.  6  and  7  show  the  same  for  a 
Wi  11 i s ton  Basin  shale. 

For  the  high  velocity  Green  River  shales  (see  Figs.  2-5),  the 
realm  of  total  reflection  (below  53°  for  the  shale  of  Figs.  2  and  3  and 
below  49°  for  that  of  Figs.  4  and  5)  remains  unchanged  by  the  addition 
of  anisotropy.  The  greatest  change  in  bottom  loss.  Figs.  2  and  4,  occurs 
in  the  range  of  angle  intermediate  between  the  critical  angles  associated 
with  the  two  wave  numbers.  The  critical  angles  are  defined  as  the  angles 
at  which  the  associated  vertical  wave  numbers  change  from  being  purely 
imaginary  to  being  purely  real,  i.e.,  the  angle  above  which  the  wave 
propagates  into  the  solid.  That  angle  associated  with  the  quasi-shear 
wave  number  q-j  is  the  first  (smallest)  critical  angle,  and  the  angle 
associated  with  the  quasi-compressional  wave  number  q^  is  the  second 
(larger)  critical  angle.  Between  the  two  wave  numbers,  anisotropic  bot¬ 
tom  loss  differs  from  isotropic  bottom  loss  by  about  0.2  dB,  and  a 
prominent  (1  dB)  lowering  of  bottom  Toss  occurs  at  the  second  critical 
angle  (about  70°),  probably  due  to  a  sampling  error  related  to 
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BOTTOM  LOSS  -  dB 


FIGURE  2 

BOTTOM  LOSS  FOR  FAST  GREEN  RIVER  SHALE 
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PHASE  OF  R12-  deg 


FIGURE  3 

PHASE  OF  REFLECTION  COEFFICIENT  FOR  FAST  GREEN  RIVER  SHALE 
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FIGURE  5 

PHASE  OF  REFLECTION  COEFFICIENT  FOR  SLOW  GREEN  RIVER  SHALE 
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1°  increments  used  in  computation.  Furthermore,  note  that  the  location 
of  this  critical  angle  is  shifted  by  almost  2° 

The  change  in  the  phase  of  the  reflection  coefficient  is  about 
2°  for  the  first  of  these  shales  (Fig.  3),  and  7°  for  the  second  shale 
(Fig.  5)  near  the  second  critical  angles,  becoming  much  less  away  from 
these  angles. 

For  the  Williston  Basin  shale  (Fig.  6),  however,  there  is  no 
range  of  total  reflection  since  the  shear  velocity  is  smaller  than  the 
water  sound  speed.  The  anisotropy  changes  bottom  loss  for  this  case  at 
angles  from  zero  to  second  critical  angle  (67°).  Of  most  interest,  how¬ 
ever,  is  a  marked  decrease  in  bottom  loss  of  about  2  dB  at  low  grazing 
angles.  The  shift  in  the  second  critical  angle  is  much  more  pronounced 
here,  approximately  8°. 

Drastic  effects  are  to  be  seen  in  the  phase  for  the  Williston 
Basin  shale  (Fig.  7).  The  phase  of  the  reflection  coefficient  for  the 
anisotropic  shale  differs  from  that  for  the  isotropic  shale  by  as  much 
as  120°  below  a  grazing  angle  of  67°.  The  two  become  practically 
indistinguishable  above  this  grazing  angle. 

In  this  section  we  have  shown  that  transverse  anisotropy  yields 
nontrivial  differences  in  bottom  loss  and  the  phase  of  the  reflection 
coefficient  when  compared  to  the  same  quantities  in  isotropic  sediments. 
The  curves  are  not  qualitatively  different,  and  most  of  the  effects  seem 
attributable  to  the  distortion  in  the  curves  due  to  a  shift  in  the  posi¬ 
tion  of  the  critical  angle  for  quasi-compressional  waves.  The  differ¬ 
ences  in  bottom  loss  amount  to  as  little  as  a  few  tenths  of  decibels 
for  low  anisotropy  sediments  (Green  River  shale),  to  as  much  as  3.5  dB 
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for  a  moderately  high  anisotropy  sediment  (Williston  shale).  Phase 
differences  range  from  2°  for  low  anisotropy  sediments  to  as  much  as 
120°  for  high  anisotropy  sediments. 

C.  Parameter  Studies  on  the  Constant 

One  problem  in  our  analysis  has  been  a  lack  of  measurements  of 
the  constant  C^3-  The  determination  of  C13  requires  the  use  of  the  dis¬ 
persion  relation  or  consistency  equation  for  the  differential  equation, 
and  a  measurement  of  the  speed  of  one  of  the  wave  solutions  to  the  dif¬ 
ferential  equations  in  some  specified  non-axial  direction.  This  non- 
axial  measurement  is  not  normally  made.  Furthermore,  a  certain  modicum 
of  error  in  C13  is  bound  to  result  in  its  computation  from  the  dispersion 
relation,  the  components  of  which  are  experimentally  determined. 

Until  measurements  become  readily  available,  it  is  of  interest 
to  consider  how  good  an  approximation  might  be  made  of  the  actual  value 
of  this  elastic  parameter.  In  this  section  we  develop  a  means  for  esti¬ 
mating  C^3  for  consolidated  (hard)  sediments  for  which  the  compressional 
wave  anisotropy  is  measured. 

A  parameterization  of  the  constant  C-^must  have  two  properties. 
It  must  reduce  to  the  isotropic  value  in  some  simple  limit,  and  it  must 
be  simpler  in  form  than  the  dispersion  relation. 

In  this  isotropic  limit  C-j3  reduces  to  x,  i.e.. 


Let  us  parameterize  C-|3  as 
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C13  =  6  ^11C33  '  2C44 


This  reduces  to  the  isotropic  limit  when  6=1  and  C11=C33 
B  is  then  defined  by  the  relation 


The  parameter 


8  = 


C13  +  2C 
/C11C33 
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If  we  divide  out  density  from  the  defining  equation  for  8,  and 
express  the  result  in  terms  of  shear  and  compressional  velocities,  we 
obtain  the  relation 


$  a 


* 


where  is  of  order  unity  and  relates  C^3  to  C^,  while  a ^  is  the  con¬ 
stant  of  proportionality  relating  C-^  to  C33. 

Experimentally  obtained  values  for  (Tables  II-IV)  show 
that  C^3  is  of  the  order  of  for  mildly  anisotropic  sediments  or  is 
some  small  multiple  of  C44  for  strongly  isotropic  sediments,  i.e.,  is 
of  magnitude  one  or  two.  The  Williston  Basin  shale  has  C-|3  approximately 
equal  to  twice  Likewise,  the  term  vt^C33  will  always  be  greater 

than  C33  for  transverse  isotropy.  A  measured  anisotropy,  of  100%,  com¬ 
puted  according  to  the  definitions  of  Chapter  II,  would  yield 
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/C11C33  =  3C33 


Therefore,  we  safely  consider  that  1«*2<3.  Thus,  the  values  of  a-j  and 
(%2 >  which  are  the  values  most  likely  to  occur  in  marine  sediments,  and 
the  fact  that  c$<Cp  demand  that  g<l. 

Let  us  then  vary  e  from  its  actual  computed  values  for  the 
sediments  available  to  us  until  we  reach  6=1,  in  order  to  determine  pre¬ 
cisely  to  what  extent  inaccuracy  in  6  affects  the  reflection  coeffi¬ 
cients.  Fortunately,  in  all  cases,  the  greatest  effect  of  variation  of 
6  occurs  in  the  vicinity  of  the  second  critical  angle.  For  the  Green 
River  shales  the  bottom  loss  curves  (see  Figs.  8  and  9)  are  undisturbed 
below  the  quasi-shear  critical  angles.  Near  the  quasi -compressional 
critical  angle  (~70°),  the  bottom  loss  curves  are  gradually  displaced 
to  a  maximum  of  about  0.3  dB  (in  Fig.  8). 

Beyond  the  quasi-shear  critical  angles,  phases  of  the 
reflection  coefficient  show  displacements  for  these  low  anisotropy  shales 
similar  to  those  they  encounter  in  a  total  reduction  to  isotropy;  about 
1°  for  the  lowest  anisotropy  shale  and  about  5°  for  the  higher  anisotropy 
shale.  We  omit  graphs  of  these  phases,  as  the  curves  so  overlay  one 
another  as  to  be  indistinguishable. 

In  the  case  of  the  Green  River  shale,  where  there  is  no  region 
of  total  reflection,  bottom  loss  at  the  lower  grazing  angle  is  affected 
less  than  0.1  dB  (see  Fig.  10),  even  though  phenomena  near  the  critical 
angle  for  quasi-compressional  waves  are  displaced  by  approximately 
0.5  dB  in  the  full  reduction  to  6=1.  However,  the  region  of  lower 
grazing  angles  where  greatest  effect  of  anisotropy  occurred  is  not 
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FIGURE  8 

GRADUAL  REDUCTION  OF  BETA  TO  1;  EFFECTS  ON 
BOTTOM  LOSS  FOR  FAST  GREEN  RIVER  SHALE 
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FIGURE  9 

GRADUAL  REDUCTION  OF  BETA  TO  1;  EFFECTS  ON 
BOTTOM  LOSS  FOR  SLOW  GREEN  RIVER  SHALE 


FIGURE  10 

GRADUAL  REDUCTION  OF  BETA  TO  1;  EFFECTS  ON 
BOTTOM  LOSS  FOR  WILLISTON  BASIN  SHALE 


ARL.UT 


changed  as  e  is  varied.  The  phases  of  the  reflection  coefficient  for 
this  shale  also  show  very  little  disturbance  as  e  is  reduced  to  one 
(Fig.  11).  The  total  disturbance  in  phase  over  the  whole  range  of  angles 
is  5°  in  the  worst  case. 

Thus  we  have  seen  that  variations  in  the  parameter  6  lead  to 

small  scale  variations  in  bottom  loss  in  the  region  of  the  critical  angle 

for  quasi -compressional  waves.  Changes  in  the  phase  of  the  reflection 
coefficient  are  also  small.  Thus  we  may  safely  make  the  approximation 
8=1,  for  consolidated  sediments.  For  hard  sediments,  therefore,  no 
great  loss  in  precision  occurs  in  using  the  approximation 

C13  =  1 1 C33  ‘  2C44 

D.  Effects  of  Anisotropy  on  Reflection  from  a  Typical  Marine  Sediment 

The  preceding  sections  dealt  with  hard  sediments  in  which  the 
densities  and  the  compressional  and  shear  velocities  were  high.  In 
typical  near-surface  silt-clay  sediment,  the  compressional  velocity  will 
be  lower  than  that  of  the  water,  and  the  shear  velocity  will  be  of  an 

order  of  magnitude  less  than  the  compressional  velocity.  It  is  of  inter¬ 

est  to  determine  whether  bottom  loss  curves  for  soft  sediments  are 
affected  in  the  same  manner  that  they  are  for  hard  sediments. 

If  we  use  a  water  velocity  of  1550  m/sec  and  a  water  density  of 

3 

1.043  g/cm  ,  we  may  use  the  ratio  of  sediment  compressional  velocity  to 
water  compressional  velocity  of  0.98  given  in  Bachman  (1979)  to  obtain  a 
sediment  compressional  velocity  of  1520  m/sec.  Use  of  the  velocity 
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FIGURE  11 

GRADUAL  REDUCTION  OF  BETA  TO  1;  EFFECTS  ON  PHASE  OF 
REFLECTION  COEFFICIENT  FOR  WILLISTON  BASIN  SHALE 


density  curve  given  in  Hamilton  (1978)  yields  a  sediment  density  of 

1.65  g/cm  .  We  use  a  shear  velocity  of  115  m/sec  (Hamilton,  1976). 

At  this  point  a  difficulty  in  interpretation  arises.  Is  the 

sediment  compressional  velocity  obtained  from  this  ratio  a  vertical 

velocity,  a  horizontal  velocity,  or  some  sort  of  mean  velocity?  Since 

Bachman  provides  no  answer  to  this  question,  we  will  compute  bottom  loss 

curves  for  each  of  the  three  possibilities:  c  =/c  c  ,  c  =c  ,  and 

P  P2  Px  P  P2 

c  =c  .  The  precise  distribution  of  velocity  between  horizontal  and 
P  Px 

vertical  is  obtained  via  the  regression  formula  of  Bachman  (1979)  for 
silt-clays,  Vy  =  0.25  +  0.83  V^,  where  V  represents  compressional  veloc¬ 
ity  in  units  of  km/sec  and  subscripts  v  and  h  denote  vertical  and  hori¬ 
zontal  directions,  respectively.  Bottom  loss  for  the  isotropic  case  is 
computed  as  a  reference.  A  summary  of  these  four  distributions  of  veloc¬ 
ity  is  contained  in  Table  V. 

With  this  information  we  can  compute  the  elastic  parameters: 

cn  =  Cp  p 

A 

r  -  .2 

C33  Pz  P 

C44  =  Cs  p 

where  p  is  the  density,  c_  is  the  horizontal  compressional  velocity, 

Px 

Cp  is  the  vertical  compressional  velocity,  and  c$  is  the  shear  velocity. 
However,  this  is  not  enough  information  to  compute  the  parameter  C-^. 
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TABLE  V 


POSSIBLE  INTERPRETATIONS  OF  CALCULATED  COMPRESSIONAL  VELOCITY 
FOR  ANISOTROPIC  AND  ISOTROPIC  SEDIMENTS 

4J 

c  c 

C  p  D 

P  KX  PZ 

_ (m/sec) _ (m/sec) _ (m/sec) 


Mean 


C  =/c  c 
p  px  pz 

1520 

1525 

1515 

Vertical 

C  “C 

p  pz 

1520 

1530 

1520 

Horizontal 

c  =  c 
p  px 

1520 

1520 

1512 

Isotropic 

1520 

1520 

1520 
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Let  us  therefore  utilize  the  parameterization  of  obtained  in  Section 
IV. C,  and  set 

C13  =  0  /C11C13  "  2C44 

For  hard  sed;  ,ients  it  was  determined  that  we  could  safely  set  3=1.  This 
proves  nothing  for  soft  sediments,  however,  since  the  magnitude  of  6  for 
soft  sediments  is  not  known.  Therefore,  bottom  loss  curves  will  be  com¬ 
puted  at  various  values  of  6  to  determine  the  effect  of  anisotropy  and 
the  sensitivity  to  0. 

In  Fig.  12  we  see  bottom  loss  computed  at  3=1  for  the  four 
distributions  of  velocity  outlined  in  Table  V.  The  principal  effect  of 
the  transverse  anisotropy  has  been  an  increase  in  bottom  loss  near  the 
angle  of  intromission  where  the  peak  in  bottom  loss  occurs.  The  actual 
value  of  bottom  loss  at  the  angle  of  intromission  is  not  significant  due 
to  increments  used  in  computation.  The  various  assignments  of  the  values 
of  compressional  velocities  lead  to  a  scatter  over  about  3°  of  the  loca¬ 
tion  of  the  angle  of  intromission.  The  assignment  cp=cp  leads  to  an 
angle  of  intromission  identical  with  that  for  the  isotropic  case.  The 

other  assignments,  c  =/c  c  and  c  =c  ,  lead  to  shifts  in  the  position 

P  Pz  Px  P  Pz 

of  angle  of  intromission  of  2°  and  3°,  respectively. 

In  Fig.  13  we  see  bottom  loss  computed  at  3=0.75  for  the  four 
distributions  of  velocity  outlined  in  Table  V.  The  cluster  of  values  of 
the  angle  of  intromission  has  displaced  about  6°  from  the  angle  of  intro¬ 
mission  for  an  isotropic  sediment.  This  displacement  of  the  angle  of 
intromission  has  so  displaced  the  curve  that  in  the  range  of  grazing  angl 


FIGURE  12 

BOTTOM  LOSS  FROM  SOFT  SILT-CLAY  SEDIMENT;  BETA  -  1 
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from  5°  to  70°  the  difference  in  bottom  loss  varies  over  the  range  from 
1  dB  (at  the  endpoints  of  the  range  of  grazing  angle)  to  55  dB  at  the 
new  angle  of  intromission,  and  is  of  the  order  of  10  dB  over  large 
reaches  of  grazing  angle. 

In  Fig.  14,  the  bottom  loss  has  been  computed  at  6=0.5  for  the 
four  distributions  of  velocity  outlined  in  Table  V.  The  angle  of 
intromission  has  shifted  even  further,  now  approximately  25°  beyond  that 
for  an  isotropic  sediment,  and  differences  in  bottom  loss  similar  to 
those  discussed  in  the  preceding  paragraph  are  now  spread  over  the  range 
of  grazing  angles  from  0°  to  80°. 

At  this  point  enough  of  a  trend  has  been  established  to  assert 
that  the  location  of  the  angle  of  intromission  is  very  sensitive  to  the 
value  of  the  parameter  g.  If,  however,  this  parameter  is  of  the  order 
obtained  for  hard  sediments,  i.e.,  somewhere  between  0.7  and  1.0,  and 
closer  to  1.0  for  less  anisotropic  sediments,  the  approximation  g=l 
obtained  in  Section  IV. C  should  hold.  The  final  verdict  on  this  point 
must  obviously  await  the  actual  measurement  of  the  elastic  parameter  C-^ 
for  soft  sediments. 
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FIGURE  14 

BOTTOM  LOSS  FROM  SOFT  SILT-CLAY  SEDIMENT;  BETA  -  0.5 
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CHAPTER  V 


CONCLUSIONS 

Previous  mathematical  formulations  of  the  problem  of  reflection 
from  the  sea  floor  have  assumed  that  the  sea  floor  is  an  isotropic 
solid.  It  was  always  known,  however,  that  this  was  at  best  an  approxima¬ 
tion.  Recent  geological  studies,  summarized  above,  have  led  to  the  view 
that  the  sea  floor  can  be  better  modeled  as  a  transversely  isotropic 
solid. 

The  equations  of  motion  for  transversely  isotropic  media 
depend  on  five  elastic  parameters  rather  than  two.  For  motion  in  only 
the  vertical  plane,  one  of  these  elastic  parameters  has  no  effect.  The 
elimination  of  this  one  parameter,  however,  also  eliminates  considera¬ 
tion  of  shear  anisotropy,  as  only  a  single  shear  velocity  occurs  in  the 
equations  of  motion  for  the  vertical  plane.  Shear  anisotropy  is  impor¬ 
tant  only  for  motion  entirely  within  a  horizontal  plane. 

The  assumption  of  plane  wave  solutions  to  the  equations  of 
motion  yields  an  equation  of  consistency  leading  to  two  solutions  for 
the  vertical  wave  number.  One  corresponds  to  quasi -compressional  motion, 
and  the  other  to  quasi -shear  motion.  An  examination  of  the  vertical  wave 
number  shows  that  it  is  possible,  in  environments  of  velocity  contrast 
sufficient  to  have  realms  of  total  reflection,  to  have  backward 
traveling  phase  fronts. 
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Using  standard  boundary  conditions,  we  solve  for  the  reflection 
coefficient  between  two  semi-infinite  media,  and  set  up  the  system  of 
equations  for  layered  media  in  a  form  amenable  to  computer  solution.  We 
show  that  the  reflection  coefficient  does  reduce,  in  the  isotropic  limit, 
to  the  usual  form. 

Studies  of  reflection  from  consolidated  anisotropic  sediments 
using  available  sets  of  parameters  were  presented.  These  show  that 
velocity  anisotropy  alters  bottom  loss  significantly,  particularly  in 
sediments  for  which  the  velocity  contrast  is  such  that  the  region  of 
total  reflection  vanishes,  or  is  greatly  curtailed.  The  degree  of  bot¬ 
tom  loss  change  from  isotropy  appears  to  be  an  increasing  function  of 
anisotropy,  and  even  for  fairly  moderate  anisotropies  is  large  enough  to 
measure.  In  highly  anisotropic  consolidated  sediments,  the  deviation 
from  isotropic  bottom  loss  amount?  to  1  dB  over  a  wide  range  of  low 
grazing  angles,  and  in  all  cases  the  position  of  the  critical  angle  for 
compressional  waves  is  shifted  more  than  1°. 

The  principal  bar  to  calculating  meaningful  results  on  the 
effect  of  anisotropy  lies  in  the  fact  that  few  full  sets  of  elastic 
parameters  have  been  measured.  We  have  shown,  however,  that  this 
deficiency  can  be  obviated  for  consolidated  sediments  by  using  measured 
compressional  wave  anisotropy  to  approximate  the  elastic  parameter  C^3 
(which  is  rarely  measured).  The  approximation  is  especially  good  if  low 
grazing  angle  phenomena  are  of  greater  concern  than  those  of  high  grazing 
angle. 

For  soft  (unconsolidated)  sediments  we  have  estimated  full  sets 
of  elastic  parameters  and  calculated  bottom  loss.  If,  as  in  the  case  of 
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hard  sediments,  a  low  anisotropy  sediment  corresponds  to  a  value  of  the 
parameter  very  close  to  its  value  in  the  isotropic  limit,  the 
approximation  discussed  above  is  valid  for  soft  sediments,  and  the  prin¬ 
cipal  effect  of  anisotropy  in  soft  (silt-clay)  sediments  is  a  small 
shift  (1-3°)  in  the  angle  of  intromission.  If  this  is  not  true,  we  have 
shown  that  bottom  loss  is  very  sensitive  to  the  value  of  3.  Obviously, 
more  information  is  needed  on  the  actual  value  of  the  constant  C-jg  in 
soft  sediments. 

In  the  work  completed  so  far,  we  have  shown  that  compressional 
anisotropy  has  significant  effects  on  both  bottom  loss  and  the  phase  of 
the  reflection  coefficient  for  both  consolidated  and  unconsolidated  sedi¬ 
ments.  The  course  of  future  work  in  this  topic  seems  obvious.  The 
assumption  of  isovelocity  layers  must  be  abandoned  and  standard  linear 
gradients  in  velocity  adopted.  This  would  include  refraction  phenomena 
known  to  be  important  in  isotropic  media.  The  computational  tools 
already  stand  in  such  form  that  the  effect  of  anisotropy  on  energy 
absorption  can  be  taken  into  account.  Additional  work  may  also  be  done 
on  a  topic  totally  excluded  from  this  paper,  that  of  surface  (interface) 


waves. 
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